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Abstract 

Metabolic Flux Analysis is now viewed as essential to elucidate the metabolic pattern of cells and to design appropriate 
genetic engineering strategies to improve strain performance and production processes. Here, we investigated carbon flux 
distribution in two Streptomyces coelicolor IKi (2) strains: the wild type IVIMS and its derivative mutant Ml 146, in which gene 
clusters encoding the four main antibiotic biosynthetic pathways were deleted. Metabolic Flux Analysis and ^^C-labeling 
allowed us to reconstruct a flux map under steady-state conditions for both strains. The mutant strain Ml 146 showed a 
higher growth rate, a higher flux through the pentose phosphate pathway and a higher flux through the anaplerotic 
phosphoeno/pyruvate carboxylase. In that strain, glucose uptake and the flux through the Krebs cycle were lower than in 
M145. The enhanced flux through the pentose phosphate pathway in M1146 is thought to generate NADPH enough to face 
higher needs for biomass biosynthesis and other processes. In both strains, the production of NADPH was higher than 
NADPH needs, suggesting a key role for nicotinamide nucleotide transhydrogenase for redox homeostasis. ATP production 
is also likely to exceed metabolic ATP needs, indicating that ATP consumption for maintenance is substantial.Our results 
further suggest a possible competition between actinorhodin and triacylglycerol biosynthetic pathways for their common 
precursor, acetyl-CoA. These findings may be instrumental in developing new strategies exploiting S. coelicolor as a platform 
for the production of bio-based products of industrial interest. 
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Introduction 

Since the late 1960s, considerable efibrt has been devoted to 
discover new antibiotics. However, the discovery of biologically 
active molecules has proved progressively more difficult due to 
biochemically and technically long and costly processes. Further- 
more, the prevalence of multi-resistant bacteria has now increased 
quite substantially [1]. Therefore, discovering novel antibiotics is 
now seen as critical for public health and medical research 
programs. 

In this regard, the soil-inhabiting. Gram-positive bacteria, 
belonging to the genus Streptomyces, are significant sources of new 
bio-active molecules [2]. While these bacteria produce up to 70% 
of antibiotics currendy used in medical prescriptions as well as 
many valuable compounds (such as immune-suppressors or anti- 
carcinogenic agents), recent publications of the genome sequence 
from several Streptomyces species [3-5] highlighted that their 
potential for the production of new secondary metabolites remains 



enormous. In fact, genome mining revealed that for a given 
Streptomyces species, the genome comprises 20 to 40 gene clusters 
presumably involved in the production of many secondary 
metabolites, while three to five only are currently produced under 
ordinary laboratory culture conditions. 

The genetic organization of genes encoding proteins involved in 
secondary metabolite production has been studied extensively. 
Several gene clusters encoding antibiotic biosynthetic pathways, as 
well as specific or pleiotropic regulators, have now been described 
[6]. 

However, although primary metabohsm provides precursors, 
reductive power and energy used for biomass and secondary 
metabolites biosyntheses, an integrated picture of relationships 
between primary and secondary metabolism is still missing. A 
metabolic flux analysis describing carbon flux distribution within 
the metabolic network, that would complete previous transcrip- 
tomic [7] and proteomic studies [8], would clarily the mechanisms 
by which primary metabolism sustains antibiotic production. This 



PLOS ONE I www.plosone.org 



1 



December 2013 | Volunne 8 | Issue 12 | e84151 



Carbon Fluxes in S. coelicolor 



would in turn be of crucial importance for microbiological 
engineering of Steptomyces strains in order to define conditions 
under which the accumulation of precursors, preferential expres- 
sion of the above-mentioned gene clusters, or antibiotic synthesis 

are optimal. 

Streptomyces coelicolor A3 (2) M 1 45 is a model strain often used to 
investigate the regulation of antibiotic production and its genome 
was the first one to be sequenced and annotated amongst 
Steptomyces species [3]. Depending on growth conditions, this 
strain may produce up to four antibiotics, but generally synthesizes 
mainly the blue pigmented, polyketide antibiotic actinorhodin as 
the major product [9]. 

Considering the production yield of S. coelicolor M145, which 
presumably indicates an efficient synthesis of the polyketide 
precursor acetyl-CoA, S. coelicolor Ml 45 has been chosen as a 
'genetic platform' to create a super-host for the heterologous 
expression of gene clusters associated with secondary metabolism 
from other Streptomyces species as well as organisms from other 
genera [10-12]. In order to a\'oid competition between native and 
heterologous biosynthetic pathways for a limited pool of metabolic 
precursors, the four major antibiotic clusters of 5'. coelicolor Ml 45 
were deleted, yielding strain Mil 46 [11]. 

In S. coelicolor M145, production of the Calcium-Dependent 
Antibiotic (CDA) coincides with growth [13], whereas production 
of actinorhodin, undecylprodigiosin and methylenomycin usually 
take place at the onset of the stationar)- phase when growth rate 
has slowed down [14,15]. In the case of actinorhodin, the shift 
between the growth and the production phases is accompanied by 
a series of transcriptomic and proteomic switches leading to 
carbon redirection to actinorhodin synthesis [7,8,16]. 

In the present paper, w(; d(;vek)ped a N-limited minimal 
medium to study S. coelicolor metabolic fluxes under steady-state 
conditions. N-limitation was chosen as a typical condition 
promoting actinorhodin production by S. coelicolor Ml 45, as 
previously described [17-19]. Under such conditions, both Ml 45 
and M 1 1 46 strains grow exponentially at a very low growth rate 
(//~0.04 h ), M145 producing only actinorhodin during its 
growth phase. 

We took advantage of metabolic steady conditions to compute a 
carbon flux distribution in S. coelicolor using a combination of 
stoichiometric inventory described by Holms [20], enzyme activity 
assays, mass-balance (Metabolic Flux Analysis) and C isotopic 
methods ('^C Metabolic Flux Analysis). Data from genome [3], 
transcriptome [7] and proteome [8] analyses were used to 
constraint and define operating metabolic pathways. 

We elucidate the flux distribution in primary carbon metabo- 
lism in both S. coelicolor M 1 45 and M 1 1 46 strains and compare 
them so as to gain insights on effects of the impairment of 
secondary metabolism on metabolic homeostasis. 



Media and growth conditions 

Batch cultures were grown in a bioreactor (Infors Labfors 3) 
with a 2 litre working volume, 2 Rushton-type impellers and 3 
baffles. The bioreactor was equipped with a cooled condenser to 
avoid medium evaporation. Temperature and pH were controlled 
at 28°C and 7.0, respectively, with automatic addition of NaOH 
(1 mol L~') or HCl (0.5 mol L~'). COg and O2 concentrations 
were monitored with a Bluesens BCpreFerm gas analyzer. 
Aeration was adjusted to 0.4 volume air (volume culture) ' 
min~' (wm) and agitation was 800 rpm so as to maintain 
dissolved oxygen near 70% saturation (aerobic conditions). 

A defined minimal medium (File SI) was developed according 
to Egli and Fiechter [21]. Based on theoretical excess factors for 
minereil elements and vitamins with respect to carbon [21] the 
medium, although N-limited, allowed S. coelicolor strains to reach a 
bacterial density of 9.6 g (cell dry wc-ight) L ' and to grow for 
more than 5 generations in the (exponential growth phase. Note 
that exponential growth during at least 5 generations was a 
prerequisite for steady-state isotope labeling in metabolites [22]. 
Isotopic labeling was conducted with a mixture consisting of 80% 
unlabelled glucose and 20% [l-'^C] glucose (CortecNet). For 
fermentation experiments, 100 mL of medium no. 1 (15.0 g L ' 
glucose, 15.0 g L ' glycerol, 15.0 g L ' bacto peptone (Difco 
0118-07), 3.0 gL"'NaCl, 1.0 gL"' CaCOs, pH 7.2) in a bafiled 
polycarbonate 500-mL Erlenmeyer flask were inoculated with 10^ 
spores from a 30% glycerol stock. The culture was incubated at 
28°C and agitated with a magnet (Nuovall agitator, Bioblock 
Scientific/Sybron thermolyne) for 48 h. Eight mL were then used 
to inoculate 100 mL of medium no. 2 (33.0 g L ' glucose, 15.0 g 
L"' yeast extract (BD 212720), pH 7.2). The second culture was 
incubated for 24 h. Cells from 40 mL were collected by 
centrifugation (3,200 g, 5 min, room temperature), suspended in 
20 mL of TES solution (270.5 mmol L~'). and immediately used 
to inoculate the bioreactor. An additional sample of 10 mL was 
used for cell dry weight quantification, after filtration through a 
cellulose acetate filter (0.22 |tm) and distilled water washing before 
drying at 95°C for 24 h. 

Gas exchange 

In bioreactors, liquid and gas phases were well mixed and gas 
diffusion was rapid such that both the CO2 Exhaust Rate [CER) 
and the O2 Uptake Rate {OUR) could be calculated at each time 
with online measurements: 

1.013*10^*60 ,„ 

CER = * f * 

103* 8.31 *(r°C-l-273.15)* 100 

0.03 



79.07 * %C02l"' 



100-%O2° 



-%C02o"' 



Materials and Methods 

Biological material 

The present study was carried out with S. coelicolor A3 (2) strain 
Ml 45 and S. coelicolor A3 (2) strain Ml 146 (act, red, cpk, cda 
deletions). S. coelicolor strain A3 (2) Ml 146 cannot produce 
actinorhodin (ACT), undecylprodigiosin (RED) and calcium 
dependent antibiotics (CDA) and no cryptic type I polyketide 
synthase (Cpk) was expressed. Both strains were a generous gift 
from M.J. Bibb and J.P. Gomez-Escribano from the John Innes 
Centre (Norwich, UK). 



OUR-- 



1.013* 10^*60 



*/«* 



103* 8.31 *(r°C-h273.15)* 100 
79.07 * %02f' 



20.9- 



100-%O2°«'-%CO2»' 
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In these equations, it is assumed that the build-up of dissolved 
CO2 and N2 in the Hquid phase was negligible. Flows ifg") as well 
as mole fractions (% CO2, % O2) were measured with dehydrated 
air to avoid any mole fraction imbalance due to water vapour. 
CER and OUR were plotted versus time and summed in order to 
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quantify total CO2 evolution and O2 consumption. The latter 
values are used in section Materials and methods 12.1. to 
determine Tco2/s and respectively. AU the abbreviations 

and symbols used in this study are listed in File S2. 

Nucleic acids, proteins, antibiotics, triacylglycerol 

measurements and elemental composition 

Nucleic acids were cjuantified after Norris and Ribbons [23] and 
proteins were quantified with the Bradford method [24]. 
Actinorhodin and undecylprodigiosin contents were determined 
spectrophotometricaUy according to Christiansen [25] and Tsao et 
al. [26] respectively. Production of CDA was detected using a 
bioassay adapted from Lautru et al. [27]. Triacylglycerol contents 
were determined by Fournier Transformed Infra Red Spectros- 
copy (FTIR) as described by Verleyen et al. [28] and Hong et al. 
[29]. FITR analyses were carried out at the Laboratoire de 
Chimie Chimie Physique, UMR 8000, CNRS/Universite Paris 
Sud, Orsay, France (see File SI for further details). The elemental 
composition (C, H, O, N) of S. coelicolor was determined with an 
elemental analyser (VarioCube, Elementar) on cells harvested 
during the exponential phase of growth. Elemental analyses were 
carried out at the Institut de Chimie des Substances NatureUes, 
UPR 2301, CNRS, Gif sur Yvette, France. 

Proteins precipitation for metabolite analyses 

Supernatants were de-proteinated by acid precipitation with 
HCIO4 35% (0.1 mL (mL supernatant) ') and then centrifuged 
2 min at 6,000 g at room temperature. Supernatants were 
neutralized with cold KOH (7 mol L~'), centrifuged 2 min at 
6,000 g at room temperature and filtered through a 0.22 \ua 
polyether sulfone membrane. 

Glucose and extracellular metabolites quantification 

Glucose, a-ketoglutarate (2-oxoglutarate), pyruvate, L-lactate, 
acetate and ethanol were analysed from de-proteinated superna- 
tants on an Agilent 1 200 HPLC instrument with an Aminex HPX- 
87H column at 35°C using H2SO4 (5 mmol L ') as the eluent 
(0.5 mL min '). Metabolites were detected either with an Agilent 
SPD UV detector (210 nm) or an Agilent RID detector (40°C). 

Ammonium quantification 

The ammonium content was measured enzymaticaUy in 20 |J,L 
of de-proteinated samples with glutamate dehydrogenase in 1 mL 
of an assay mixture made of Tris/HCl (0.1 mol L~') pH 8.0, 
NADH (0.2 mmol L"'), ADP (0.6 mmol L"'), a-ketoglutarate 
(10 mmol L~'), and 5 units of glutamate dehydrogenase from 
bovine liver (Sigma #G2626). The ^340 was monitored until the 
end of the reaction (~8 min). 

Phosphate quantification 

Phosphate was measured on 50 |J,L from de-proteinated samples 
using 100 |xL of PiBlue (PiBlue Phosphate Assay Kit, Centaur). 
The ^600 was monitored for 30 min at room temperature. 

Enzymatic assays 

40 mL of culture were harvested and bacterial cells were 
suspended in 10 mL Tris-Tricarballylic acid (15 mmol L ') buffer 
pH 7.8, MgCla (10 mmol L"'), 10% glycerol. RNase A (1 mg 
mL ') and DNase I (1 mg mL ') were added when needed. The 
sample was sonicated 8 times with a Branson SoniEer 250 
sonicator set at a power level = 2, during 20 seconds at 30 seconds 
intervals. Cell debris was removed by centrifuging at 10,000 g for 
15 min at 4°C. The supernatant was used as the crude extract. 



The protein concentrations were measured by the Bradford 
method using BSA as a standard [24] . All enzyme activities were 
measured on a Beckmann DU7400 spectrophotometer at 28°C 
according to Sugimoto and Shiio [30,31], Nachlas et al. [32] and 
Sauer et al. [33] with slight modifications (see File SI for further 
details). 

Sample preparation and GC-MS analyses 

GC-MS analyses of amino acids to determine the ''^C-labeling 
(see File SI for further details) were carried out as described by 
Zamboni et al. [34]. Briefly, after 5 generations of exponential 
growth in the presence of 20% [l-'^^C] glucose, cells were 
collected, proteins were extracted and submitted to acidic 
hydrolysis, extracts were derivatized with methoxamine and N- 
methyl-N-(trimethylsilyl) trifluoroacetamide and injected on a GC- 
MS (Agilent 689()N gas chromatograph, time-of-flight mass 
spectrometer Pegasus III). Peak integration was performed using 
LECO Pegasus software. Because automated peak integrations 
have been occasionally unreliable, manual controls and/ or 
corrections were systematically performed for each metabolite. 

Five independent samples were collected along a 1.5 -generation 
period and analysed in triplicate. During the sampling period, ji, 
qs, qo2, yco2, ^x/s, I'act/s, Yco2/s and ¥02/8 values were constant 
and at their maximal vEilues. MID values for Glu, Thr, Asp, Ala 
and Val were identical, with less than 1% RDS. We therefore 
conclude that both metabolic and the isotopic steady-states were 
ensured in our experiment. 

S. coelicolor central metabolism 

Pathways considered. The central metabolic network of S. 
coelicolor was deduced primarily from genome annotation [3], 
transcriptome [7] and proteome analyses [8] . 

Glucose c-ntry and phosphorylation through the Glk-GlcP 
permease complex in S. coelicolor were well described [35—37]. 
Consequendy, we assumed that glucose phosphorylation was 
carried out by the ATP-dependent glucose kinase Glk. 

All genes and enzymes involved in the Embden-Meyerhof- 
Parnas pathway (EMP), the tricarboxylic acid cycle (TCA), and the 
pentose phosphate pathway (PPP) were detected in S. coelicolor 
[8,38-40]. 

The Entner-Doudoroff pathway (ED) has been shown to be 
non-functional in several Streptomyces species [39,41]. Since (i) there 
was Utde growth with two carbon substrates and (ii) both isocitrate 
lyase and malate synthase were not detected in the S. coelieohr 
proteome [8], the glyoxylate shunt was considered to be negligible 
in our model. 

In the present paper, comparison of the experimental and 
calculated yco2 values, and analysis of '"^C-labeling in valine, 
alanine, glutamate, glutamine and proline provided further 
evidence of these assumptions for S. coelicolor under our experi- 
mental conditions (see 3.4.). 

PEP carboxylase was shown to be the sole anaplerotic enzyme 
in Streptomyces species [8,38,39,42,43] and was therefore considered 
in our model. 

Few experimental data on the respiratory chain efficiency in 
Streptomyces species are available [44]. From consideration of 
genome mining, we assumed here that the respirator^' chain of S. 
coelicolor could be branched, with a type I NADH dehydrogenase 
excreting 4 protons per NADH, an uncoupled type II NADH 
dehydrogenase and two terminal oxidases, the bd complex that 
excretes 2 protons per NADH and the bci/aag supercomplex that 
excretes 6 protons per NADH [45]. The t)'pe I NADH 
dehydrogenase has been shown to account for about 70% of 
NADH dehydrogenase activity in aerobic conditions in Escherichia 
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coli [46]. The supercomplex bci/aas was shown to have a lower 
afiinity for oxygen and to be the predominant complex in aerobic 
conditions [47]. Enzymes corresponding to both type I NADH 
dehydrogenase and to hcj/aa:^ supercomplex were identified in S. 
coelicolor proteome [8] . Genome mining and proteome analyses [8] 
further revealed the presence of an FiFq-ATP synthase. Under the 
assumption that the synthesis of one ATP by FiFq-ATP synthase 
requires three to four protons, NADH re-oxidation via the bcl/aaS 
branch yields 1 .50-3.33 ATP. FADH re-oxidation was assumed to 
be due to succinate dehydrogenase [8,44]. Depending on the 
pathway used by electrons, the transfer of two electrons from 
succinate to oxygen yields 0.50-2.00 ATP [48]. 

Considering all these data, a P/O value (moles of ATP 
produced per 0.5 mole of O2 consumed) of 2.0 could be 
reasonably considered for S. coelicolor. Nevertheless, our results 
on carbon flux distribution are not sensitive to the P/ O value used. 

Two genes of S. coelicolor {SC07622 and SC07623) putatively 
encode a PntAB-type nicotinamide nucleotide transhydrogenase 
that interconverts NADP + NADH to NADPH + NAD in E. coli 
[33], and thus the presence of a transhydrogenase was taken into 
account. 

Amino acids, nucleic acids, peptidoglycan, lipids and 
small molecules biosyntheses. Amino acids, nucleic acids, 
small molecules, peptidoglycan, triacylglycerol and lipids biosyn- 
thetic pathways in S. coelicolor appeared to be identical to those 
present in other bacteria [42] and therefore their stoichiometry 
was assumed to be the same as that observed in other 
microorganisms [22,49-61] (see File S3). 

Antibiotic biosyntheses. The biosynthetic pathways for the 
3 known antibiotics produced by S. coelicolor were reconstructed 
using S. coelicolor genome and published data [62-70]. For each 
antibiotic, all the steps of the biosynthesis were fused into one 
single reaction that includes precursors, reducing power and 
energy needs. 

Biomass composition. To cak:ulate carbon fluxes according 
to the stoichiometric approach of Hokns [20], the biomass 
composition of S. coelicolor cells had to be determined. The 
respective amounts of the main polymers involved in the S. 
coelicolor biomass composition, namely proteins, DNA, RNAs, 
triacylglycerol, lipids and peptidoglycan were taken from Borodina 
et al. [71] with some modifications according to the experimental 
values we determined. These modifications are discussed in section 
Results 1. 

Precursors, cofactors and energy requirements. Using 

the biomass composition, reactions assot:iated with building-block 
biosynthesis (primary metabolism) and reactions for secondary 
metabolism, we calculated precursor, cofartor and energy 
requirements to synthesise 1 g of dry cells. On the basis of 
Ingraham's calculations for E. coli [22], the energy requirement for 
the polymerisation of macromolecules in S. coelicolor was calculated. 
We also calculated specific precursor, energy and reducing power 
requirements to synthesize 1 mole of antibiotics. 

Flux calculations. 

Experimental kinetic parameter quantifications. For 

each experiment, we calculated the following graphs: ln(X(t)) 
versus t, X(t) versus S(t), ACT(t) versus S(t), RED(t) versus S(t), 
KETft) versus S(t), PTRft) versus Sft), COzftJ versus S(t), O^ft) versus 
S(t), where X denotes biomass and S denotes source glucose. 
During exponential growth, each curve was linear, and the slope 
allowed us to determine growth rate constant n (h~') and the 
following yields: Tx/s, Yact/s, Y red/s, 'i'cDA/s, I^kkt/s, Ynn/s, 
^C02/s, ^02 /s- For example, Tx/s is the yield of biomass production 
with respect to glucose, etc. Note that yields are expressed in 



carbon equivalents, e.g. Tx/s is in moles of biomass carbon per 
mole of glucose carbon. Using these values, the glucose 
consumption and productivity of metabolites (jlmol (g cell dry 
weight) ' h ') were computed as (note that the following are 
algebraic and therefore can be negative): 

Yx/s*MWx*6 



_ Y ACT IS * Is 
7 ACT — ^ 



_ Yket/s * Qs 
Vket — ? 



_ YpYR * qs 

YPYR — 5 



'yC02= - {Yc02/S * qs) 



qm = Yoi/s * qs 



Carbon balance (in %) was calculated as follows: 



^balance 

(Xx/s +Yco2is +Yactis +Yred/s +Ycda/s+Ypyr/s+Yket/s) * 10^ 

Anabolic fluxes calculation. Glucose-6-phosphate, fructose- 
6-phosphate, ribose-5-phosphate, erythrose-4-phosphate, glyceral- 
dehyde-3-phosphate, 3-phosphoglycerate, phosphomoi^yruvate, 
pyruvate, acetyl-CoA, a-ketoglutarate and oxaloacetate are the 
11 precursors from central metabolism. These precursors are 
sufficient to synthesize any macromolecule in the bacteria. That is, 
carbon flux to anabolic reactions can be represented by 1 1 
individual fluxes from precursors to biomass and/or secreted 
metabolites (antibiotics, organic acids). Known molar require- 
ments of biomass or antibiotic for precursors (for a given precursor 
i, they are denoted as pi for biomass, /, for actinorhodin, m,- for 
undecylprodigiosin, n, for the calcium dependent antibiotic 
expressed in |imol (g ceU dry weight) and |imol (|imol 
antibiotic)"', respectively) and productivity {q and y values, see 
Materials and Methods 12.1) allow one to calculate anabolic 
requirement fluxes i>i and the sum ii,- for any precursor i, as 
follows: 

'Pf =Pi* li (biomass) 



(pf^^ = yact * h (actinorhodin) 



0RED _ ^ ffj. (undecylprodigiosin) 
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=ycDA* (calcium) 



-Vi (organic) 

^NA = ^jS- + 0ACT _^ ^RED _^ ^CDA _^ ^SEC 

Metabolic Flux Analysis. Although anabolic fluxes were 
calculated as described in Materials and Methods 12.2, catabohc 
fluxes in central metabolism CEilculations remained to be deter- 
mined to obtain a full carbon-flux distribution. For this purpose, 
metabolic reactions were represented as a stoichiometric matrix [S) 
of size {m, n}. Each row of this matrix represented one unique 
compound amongst the m = 48 compounds and each column 
represented one reaction amongst the n = 46 lumped metaboUc 
reactions. Note that the 46 lumped reactions correspond to 265 
enzymatic reactions. S was a sparse matrix, with 0 as the most 
represented value. The flux through all reactions in the metabolic 
network was represented by the vector v (dimension n). The set of 
concentrations in metabolites was represented by the vector x 
(length m). The system dynamics is described by dx/dt = Sv. In the 
steady state, dx/dt=0, that is, Sv = 0. In the present study, this 
equation is true as metabolism was studied at the steady state, that 
is, in the exp<)n(;ntial growth phas(' in which turn-over fluxes are 
constant and therefore there is no build-up in intracellular 
metabolite concentration. Glucose specific consumption (qs) and 
the 1 1 anabolic fluxes (CP;) are used as external constraints. 
However, the system stiU has one degree of freedom. Data from 
labeling experiments were used to solve the system. O2 consump- 
tion, CO2 evolution and the respiratory quotient [QR) were not 
used as constraint here and were rather compared to the 
computed O2 and CO2 fluxes to check the reliability of our 
system, as discussed in section Results 2. 

'^C Metabolic Flux Analysis. Here, we took advantage of 
the software OpenFlux, which is a computational tool based on 
the Elementary Metabolite Units (EMU) framework [72,73] to 
analyse '^C-labeling data. The EMU framework used here for the 
simulation is provided in Table SI. OpenFlux calculates fluxes 
from Mass Isotopic Distribution (MID, or isotopic pattern) of 
amino acids. [''*C-1] glucose was used here as the labeling 
substrate and allowed us to determine the flux-ratio between the 
pentose phosphate pathway and glycolysis (PPP/Glyc ratio). In 
fact, when ['''C-1] glucose was metabolized by the PPP, the '^C 
label was lost by decarboxylation, while the ''^C label was 
incorporated into amino acids when it enters glycolysis. We 
recognize that ri'\'ersible triose phosphati; isomerase and aldolase 
activities could lead to a '''C-labeUng in the C-6 atom position in 
fructose- 1, 6-bisphosphate and then possibly in hexose phosphates 
[74]. If significant, this effect should be taken into account to 
calculate the relative commitment of PPP and glycolysis. The ' C- 
labehng in tryptophan and histidine should provide evidence of the 
labeling in C-6 because their -COOH group comes from this very 
atom position. While tryptophan was absent from the spectra since 
it is degraded to a considerable extent by acid hydrolysis [75], the 
'"'C-labeling pattern in histidine (Table S2) showed that the C-6/ 
C-1 isotopic exchange is negligible. Thus, the output of our 
calculations represent actual net fluxes. 

MID data from mass-spectrometer were normalized to total 
mass-signals as follows, giving relative isotopic abundance: 



^normalized 0 

If) — 

mo+mi+m2 



^normalized _ 

' mo+ini+m2 



^normalized ^2 

mo+mi +m2 

Normahzed MID were also corrected for natural isotope (C, H, 
N, O, Si) abundance during flux calculations. The glucose specific 
consumption qs, the 1 1 anabolic fluxes 0, and the normalized 
MID in glutamate (mo/ z = 348), threonine {nio/ z = 320), aspartic 
acid (mo/ Z = 334), alanine (mo/^ = 218) and valine (in^)/ z = 24:6) 
were combined to provide a unique solution for all carbon fluxes, 
thereby solving vector v. 

For sensitivity analyses, carbon fluxes were estimated by 
minimizing the variance-weighted sum of squared residuals 
(SSR) between the simulated and model-predicted MIDs using 
non-hnear least-squares regression [76]. In all cases, flux 
estimation was repeated at least 20 times starting with random 
initial values for all fluxes to find a global solution. At convergence, 
standard deviations and 95% confidence intervals for all fluxes 
were calculated using the parameter continuation technique 
described by Antoniewicz et al. [76]. This technique is based on 
evaluating the profile of SSR as a function of one flux, while the 
values for the remaining fluxes are optimized. The 95% 
confidence interval for an evaluated flux corresponds to flux 
values for which SSR increased by less than 3.84 [76]. A flux was 
considered non-observable if the 95% confidence interval was 
larger than the estimated flux value. 

Metabolic fluxes were expressed as absolute (jlmol (g cell dry 
weight)"' h~') or relative fluxes with respect to glucose uptake 
(|imol (|imol glucose) '). Relative fluxes allowed the c()m])aris()n 
between strains, as Ml 45 and Ml 146 did have different glucose 
uptake rates. 

Energy, reducing power and gas exchange fluxes 
calculations from carbon-flux distribution. With the for- 
mula used for precursors, the stoichiometry of reactions, including 
those for antibiotics syntheses (File S3), and the flux values 
(Materials and Methods 12.4), it was possible to calculate 'intrinsic 
energy production' (ATP), reducing power (NADH, NADPH) and 
gas exchange (O2, CO2). The total CO2 evolution rate was 
computed from catabolic CO2 production (glycolysis, TCA, PPP) 
and anabolic production as: 

.MODEL _ ,^NA , „CAT 

'cm — ^cm. ^co2 

Then, this value was compared with our experimental CO2 
production rate to verify the validity of the carbon-flux distribution 
(see section Results 4). In other words, experimental CO2 
evolution was an independent data point not used as a constraint 
to calculate carbon-flux distribution. We assumed that NADPH 
requirement and production should have been equal since the 
metabolism was in steady state. If the production of NADPH was 
more important than the need, excess NADPH might have been 
converted to NADP^, thereby reducing NAD"^ into NADH, by 
nicotinamide nucleotide transhydrogenase. 
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Table 1. Biomass composition of 5. coelicolor Ml 45 cells {|imol (g dry cells) 





PROTEINS 




PHOSPHOLIPIDS 




SOLUBLE MOLECULES 




Alanine 


556 


Phosphatidylethanolamine 


29 


NAD 


6 


Arginine 


145 


Phosphoglycerol 


6 


NADP 


5 


Asparagine 


162 


Cardiolipin 


1 


CoA 


5 


Aspartate 


160 


C14 


10 


Menaquinone 8 


4 


Cysteine 


68 


C15 


31 


Tetrahydrofolate 


8 


Glutamate 


157 


C16 


28 


FMN 


8 


Glutamine 


156 


C17 


31 


FAD 


5 


Glycine 


830 


C18:l 


6 


ATP 


1.5 


Histidine 


54 


TRIACYLGLYCEROLS 




ADR 


0.5 


Isoleucine 


196 


Glycerol-3-Phosphate 


22 


PEPTIDOGLYCAN 




Leucine 


305 


C14 


1 


udp-nam" 


no 


Lysine 


205 


C15 


30 


UDP-NAG" 


132 


Methionine 


91 


C16 


7 


Alanine 


209 


Phenylalanine 


98 


C17 


28 


Diaminoplmelic acid 


125 


Proline 


'^^^^ 173 


C18:l 


0.2 


Glutamate 


11 


Serine 


189 


TEICHOIC ACIDS 




Glycine 


107 


Threonine 


191 


Teichoic acid 


16 


Bactoprenol Phosphate 


0.3 


Tryptophan 


24 


Lysine 


4 


DNA 




Tyrosine 


65 


RNA 




dAMP 


27 


Valine 


328 


AMP 


100 


dCMP 


70 


CELL WALL 




CMP 


172 


dTMP 


27 


UDP-NAG" 


83 


TMP 


110 


dGMP 


70 


Galactose 


167 


GMP 


138 







"UDP-N-Acetylglucosamine 

""UDP-N-Acetylmuramic acid 

doi:l 0.1 371/journal.pone.00841 51 .tool 



Similarly, the metabolic requirement for NADH should match 
NADH production; but here and quite generally, NADH 
production was larger than NADH oxidation. In order to satisfy 
our steady state condition, we assumed that excess NADH was re- 
oxidized by the respiratory chain. For O2 consumption, tlie 
anabolic demand was completed by the extra consumption 
required for excess NADH re-oxidation with the following 
stoichiometry: 1 mole of re-oxidized NADH consumed 0.5 mole 
of O2. For ATP as well, anabolic and catabolic fluxes were 
completed by the production of ATP associated with excess 
NADH re-oxidation. At this stage, knowledge of the actual P/ O 
value (moles of ATP generated per mole of oxidized NADH) 
would be ideal. As stated in section Materials and Methods 11.1, 
we assumed that under 70% dissolved O2, the bcj/aay cytochrome 
pathway probably prevailed, and so a P/O value of 2 could be 
used to estimate ATP production from oxidative phosphorylation. 
Under this assumption, as for NADH, ATP production would be 
higher than ATP consumption. Excess ATP production would be 
likely compensated for by the consumption required for cell 
maintenance. Such a consumption is hereafter denoted as mATP, 
that is, energy requirements for homeostasis and turn-over of 
structural molecules as well as others processes that could not be 
quantified precisely {e.g. energy needed to excrete antibiotic). 



Results 

Composition of S. coelicolor biomass 

As mentioned in section Materials and Methods 1 1 .4, we 
mainly used here metabolite contents reported by Borodina et al. 
[71] with significant modifications for deoxynucleotides (for 
reasons that wiU become apparent below), the inclusion of 
bactoprenol (0.34 |j,mol (g cell dry weight)" ) in cell wall 
composition [50], and a modification of the triacylglycerol content 
in Ml 146. The biomass decomposition in building blocks is shown 
in Table 1 . Using the elemental formula of individual metabolites 
and taking into account the loss of water during polymerization of 
macromolecules, the building block composition allowed us to 
recalculate the macromolecular composition of bacteria in 
proteins, DNA, RNAs and lipids (Table 2). Our estimates were 
in agreement with published values [77-80] except for phospho- 
lipids and DNA. 

The dilference in phospholipid content was considered to be 
acceptable (0.42% of total biomass) since it was previously shown 
that such a low difference in biomass composition had a negligible 
impact on carbon flux determination [81]. 

The experimental DNA content obtained here was 6.0% 
(0.06 g (g cell dry weight) '). As the mass of the S. coelicolor genome 
is 1.86-fold that oi E. coli, which has a well-accepted DNA content 
of 3.1% [82], the value of 6% is much more accurate than that 
previously pubhshed for S. coelicolor [79] . 
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Table 2. Macromolecules content of S. coelicolor Ml 45 cells. 



Macromolecules 


Published data^ 


Calculated data** 


RSD"= 


Proteins 


41.2 


41.1 


0.17 


DNA 


3.6 


6.0 


35.36 


RNA 


16.7 


16.7 


0.00 


Phospholipids 


2.7 


3.3 


14.14 


Triacylglycerols 


1.8 


2.0 


7.44 


Soluble molecules 


3.0 


2.6 


10.10 


Peptidoglycan 


11.0 


11.4 


253 


Cell Wall 


4.4 


4.4 


0.00 


Teichoic acids 


6.6 


6.8 


2.11 


Ash 


9.0 


5.7 


31.75 


Total 


100.0 


100.0 


0.00 



^Biomass composition (% dry cell weight) from [71]. 
"^Calculated from our composition of building blocks (% dry cell weight) 
'^Relative standard deviation {%) between published and calculated values. 
doi:l 0.1 371 /journal.pone.00841 51 .t002 



The chemical elemental formula calculated from the data in 
Table 1 (with a correction for water content from 4 to 5% H2O 
excess [83]) was GHi 52 O0.45 Nq 25 (normalized with respect to 
carbon), with 5.7% of ash. This is close to those found in other 
microorganisms [78,84—87] and in good agreement with the 
elemental formula determined experimentally by combustion/ 
reduction-based elemental analysis. We thus considered that the 
biomass composition described in Table 1 was accurate. 

A difference in the triacylglycerol content was found in the 
biomass composition of M 1 45 and M 1 1 46, with a triacylglycerol 
content of 0.02 g (g cell dry wei ght)"' for M145 and 0.03 g (g cell 
dry weight)"' for Ml 146 (M. Vergne, personal communication). 

Using (z) the biomass composition (Table 1, Table 2), [iij the 
stoichiometry of reactions describing the biosynthesis of metabo- 
lites (File S3) and (m) the energy requirements re-calculated for S. 
coelicolor on the basis of published energy requirements for the 
polymerisation of macromolecules in E. coli [22], we calculated 
precursors, cofactors and energy requirements to synthesise 1 g of 
dry cell matter (Table 3). 

Enzymatic activities in S. coelicolor 

In Streptomjces lividans and other Actinomycetes, glucose-6- 
phosphate dehydrogenase was shown to be NADP-dependent 
whereas 6-phosphogluconate dehydrogenase was shown to be 
NAD-dependent [39,88]. Because of its importance for redox 
mass-balance, the nature of the cofactors was verified here. We 
found that glucose-6-phosphate dehydrogenase was strictiy 
NADP-dependent, with a specific activity of 144.6±2.7 |j.mol 
min ' (mg protein) ' and no detectable activity with NAD. 
Furthermore, we found that 6-phosphogluconate dehydrogenase 
was strictly NADP-dependent, with a specific activity of 
236. 4± 6.4 ixmol min ' (mg protein) ' and no detectable activity 
with NAD. We also verified that isocitrate dehydrogenase was 
strictiy NADP-dependent, with a specific activity of 
132.7±4.3 |J.mol min~' (mg protein)"' and no activity at all with 
NAD. 

The S. coelicolor genome contains two genes {SC07622 and 
SC07623) encoding a PntAB-like nicotinamide nucleotide traiis- 
hydrogenase, but no UdhA-eiicoding gene. In E. coli, PntAB 
reduces NADH using NADP as a reductant in vivo while UdliA 
catalyses the reverse reaction [33]. By enzymatic measurement, we 



Table 3. Precursors, energy and reducing power needs for 1 
g dry biomass syntliesis in S. coelicolor. 





Precursors, Energy, 
Reducing Power 


^mol (g 


cell dry mass) ^ 




M145 


IV11146 


Glucose-6-phosphate 


327 


327 


Fructose-6-phosphate 


357 


357 


Ribose-5-phosphate 


1,344 


1,344 


Erythrose-4-phosphate 


199 


199 


Glyceraldehyde-3-phosphate 


281 


292 


3-Phosphoglycerate 


1,701 


1,701 


Phospheno/pyruvate 


508 


508 


Pyruvate 


2,600 


2,629 


Acetyl Coenzyme A 


2,043 


2,234 


of-Ketoglutarate 


754 


754 


Oxaloacetate 


2,773 


2,802 


CO2 


-1,433 


-1,491 


O2 


17 


17 


ATP 


34,606 


34,849 


NADPH 


13,491 


1 3,993 


NADH 


-3,387 


-3,415 



Positive values correspond to metabolite consumptions and negative values 
correspond to metabolite productions. 
doi:1 0.1 371 /journal.pone.00841 51 .t003 



found a clear transhydrogenase activity, that catalysed the 
reduction of thio-NAD using NADPH with a specific activity of 
90.4±3.7 |a.mol min ' (mg protein) ', and no activity at all with 
thio-NADP and NADH. Since no evident NADPH phosphatase- 
encoding gene was found in the S. coelicolor genome, transhy- 
drogenase was then introduced in our model as the enzyme that 
balanced NADPH overproduction. 

Growth and metabolite excretion of S. coelicolor strins 

Batch-cultivation of the strains M 1 45 and M 1 1 46 yielded a 
bacterial density of 8 g (cell dry weight) L ' and 9 g (cell dry 
weight) L~', respectively. The stationary growth phase began 
when the medium turned to be nitrogen-depleted, glucose and 
phosphate being limiting approximately 9 h after nitrogen 
depletion for both strains. 

During the exponential phase in bioreactors, the growth rate 
constant (ji) of Ml 146 was 13 ±0.1% higher compared to Ml 45 
while the glucose specific uptake rate was nearly 13±0.5% lower 
(Table 4). The respiratory activity (O2 consumption and CO2 
production) was 20 ± 2 % lower in M 1 1 46 compared to M 1 45 and 
furthermore, strain M145 produced 2.03±0.08 [imol (g cell dry 
weight) ' h ' actinorhodin (up to 400 mg L ' at the end of the 
culture) while no antibiotic production was detected at all in 
Ml 146. 

The carbon balance during the exponential phase was similar in 
both strains (92±2% in M145 and 87± 1% in Ml 146) and slightiy 
lower than 100%, suggesting that a certain carbon amount was not 
recovered. HPLC, LC-MS analyses and bioassays on culture 
supernatants did not allow the identification of extracellular 
metabolites that could account for this apparent loss of carbon. 
Nevertheless, these carbon balances are similar to the 94±9% 
value found by Borodina et al. [41] for S. coelicolor Ml'i5 and to the 
90% value found by Bruheim et al. [89] for S. lividan.s during 
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Table 4. Kinetic pa 


rameters and GC-MS 


analyses. 








5. coe/ico/orN\^A5 


5. coelicolor M11 46 




Values 


RSD (%) 


Values 


RSD (%| 




0.0394 


0.54 


0.0446 


0.95 




qs ^imol (g cell dry 
mass)"' h"' 


664 


3.24 


575 


1.40 


qo2 fitnol (g cell dry 
mass)"' h"' 


2,235 


2.97 


1,528 


6.67 


yco2 limol (g cell dry 
mass)"' h"' 


-2,295 


3.24 


-1,372 


0.86 


•McT nmol (g cell dry 
mass)"' h"' 


-2 


5.87 


0 


0.00 




Glu-TMS (m/z= 348) mo 


0.54614 


0.24 


0.57610 


0.30 


Glu-TMS (m/z= 349) m, 


0.30993 


0.79 


0.28750 


0.34 


Glu-TMS (m/z= 350) mj 


0.14393 


0.97 


0.13640 


0.83 


Thr-TMS {m/z = 320) mo 


0.57694 


0.11 


0.60627 


0.14 


Thr-TMS (m/z = 321) m. 


0.29529 


0.31 


0.27454 


0.45 


Thr-TMS {m/z = 322) mj 


0.12776 


0.53 


0.11919 


0.98 


Asp-TMS (m/z = 334) mo 


0.57730 


0.22 


0.60477 


0.08 


Asp-TMS (m/z = 335) m. 


0.29536 


0.22 


0.27410 


0.37 


Asp-TMS {m/z = 336) mj 


0.12733 


0.52 


0.12113 


0.96 


Ala-TMS (m/z = 21 8) mo 


0.71483 


0.19 


0.72507 


0.15 


Ala-TMS (m/z = 21 9) m, 


0.20701 


0.51 


0.19744 


0.52 


Ala-TMS (m/z = 220) mj 


0.07817 


1.72 


0.07749 


0.85 


Val-TMS {m/z = 246) mo 


0.64474 


0.36 


0.66780 


0.25 


Val-TMS (m/z = 247) m. 


0.26438 


0.74 


0.24531 


0.70 


Val-TMS (m/z = 248) mj 


0.09088 


0.45 


0.08689 


0.62 



Kinetic parameters were experimental values obtained as described in 2.19.1. 
Mass Isotopic Distributions were normalized as described in section 2.19.4. 
Values were means of two independent experiments and RSD their relative 
standard deviations. Glu: Glutamic acid; Thr: Threonine; Asp; Aspartic acid; Ala; 
Alanine; Val; Valine. 

For qs, Yact, qo2 3nd yco2< positive values correspond to metabolite 
consumption and negative values correspond to metabolite production. 
doi:10.1371/journal.pone.0084151.t004 

actinorhodin production in N-limited conditions. The carbon 
imbalance probably indicates the synthesis of one or several yet 
unidentified secondary metabolites in S. coelicolor under N- 
limitation. 

During the stationary phase, both strains produced ot-ketoglu- 
tarate (82 |xmol±4 (g cell dry weight)"' h~' in M145 and 
162 |imol±6 (g cell dry weight) ' h 'in Ml 146) and pyruvate 
(21 nmol±l (gceU dry weight)"' h"' inM145 and 167 nmol±ll 
(g cell dry weight)"' h"' in Ml 146). There was, therefore, a clear 
difference in carbon allocation between the two strains, with 
M 1 1 46 liberating more pyruvate and a-ketoglutarate in the 
medium, likely as a result of carbon re-direction in the absence of 
growth and antibiotic production, as discussed in section 4.4. 

Flux patterns of S. coelicolor strains 

The use of stoichiometric coefficients and biomass composition 
to infer metabolic fluxes (Metabolic Flux Analysis) could not be 
carried out properly since an uncertain flux value remained at the 
PPP/ glycolysis branching-point. We thus used C-labeling ([ ''C- 
1] glucose) and '''C Metabolic Flux Analysis in the isotopic steady- 



state to estimate the relative commitment to PPP as opposed to 
glycolysis. Isotopic patterns obtained by GC-MS in amino acids 
from protein hydrolysates are shown in Table 4. In the amino 
acids considered (Glu, Thr, Asp, Ala and Val), there was a lower 
'■'C-labeling (mass nio proportionally more abundant) in the strain 
Ml 146, clearly showing a lower commitment of the ' 'C-label to 
glycolysis and the TCA (from which those amino acids are 
synthesized) and thus a larger commitment to PPP. 

Analysis of ''^C-label patterns in alanine and valine clearly 
showed the absence of ' 'C label in the C-1 atom position in both 
amino acids (Table S2). This demonstrated that the C-1 atom 
position in pyruvate was not labelled and thus the absence of an 
active Entner-Doudoroff pathway under our conditions. More- 
over, a discrepancy between yco2 and y^oi^^ would have been 
detected if either the Entner-Doudoroff pathway or the glyoxylic 
shunt had been active in our conditions, but this is clearly not the 
case (see also below). Similarly, the '^C-labeling in C-1 in 
glutamate, glutamine and proline showed the activity of the 
anaplerotic phosphoeno/pyruvate carboxylase and examination of 
the ' ^C-labeling at the C-5 position in glutamate, glutamine and 
proline indicated that the glyoxylic shunt was inactive (Table S2). 

Using further the data of Table 4, fluxes were numerically 
estimated with OpenFlux, which is based on optimization of 
weighted least squares (see section Materials and Methods 12.4). 

We analyzed the statistical significance and confidence of 
estimated fluxes in S. coelicolor M145 and Ml 146, with [''^C-1] 
glucose as a tracer and the sensitivity test described by 
Antoniewicz el al. [76]. Small sensitivities, i.e., large changes in 
the flux value resulting in small changes in the minimized sum of 
squared residuals, indicate that the flux cannot be estimated 
precisely. Large sensitivities, on the other hand, indicate that the 
flux is well determined. In the present study, sensivity tests fulfilled 
large sensitivity (96%) with a 95% confidence level in all cases. 
According to Wiechert et al. (1997), our flux maps are represen- 
tative descriptions of the metabolic processes in S. coelicolor [90] . 

In Ml 45 and Ml 146, the flux through the PPP was found to be 
22±0.3% and 38±0.2% of available glucose-6-phosphate, 
respectively. 

Our results from Metabolic Flux Analysis and isotopic labeling 
were combined to draw a carbon flux map in M 1 45 and M 1 1 46 
during the exponential phase of growth (Fig. 1, Fig. 2). Fluxes were 
expressed as |imol (g cell dry weight) ' h ' or normalized to 
glucose uptake (|imol (|imol glucose) '). Consistently, the com- 
mitted step of PPP (glucose-6-phosphate dehydrogenase) was 
found to be associated with a flux of 0.372±0.002 |J.mol (|imol 
glucose)"' in Ml 146, which is 1.7-fold higher than in Ml 45 
(0.214±0.003 |J,mol (|imol glucose) '). In contrast, the commit- 
ment to the TCA was 1.019±0.02 jxmol (|imol glucose)"' in 
M145, that is, 1.4-fold higher than in Ml 146 (0.730±0.018 |imol 
(|imol glucose) '). Precursor abstraction from primary carbon 
metabolism to biomass synthesis was larger in M 1 1 46 (total 
1.019±0.018 nmol (|imol glucose)"', giving 60±1% (molC 
mol 'C) of carbon uptake allocated to biosyntheses) compared 
to M145 (0.813±0.02 i^mol (nmol glucose)"', 47±r/o). There 
was a much larger commitment of phosphoozo/pyruvate to 
carboxylation by the PEP carboxylase in M 1 1 46 and therefore, 
a larger anaplerotic oxaloacetate production to sustain the TCA. 

Summing anabolic and catabolic fluxes for COj gave the total 
calculated CO2 evolution rate [Jcm'^^) (Fig- 1, Fig. 2, right hand 
side). Our calculated fluxes may be compared to experimental 
values (Table 4, last line). There was good agreement between 
them, with a diflerence of 3.8±1% (Ml 146) and 2±1% (M145) 
only. Similarly, the total O2 consumption rate {loi"^^^) '^^s 
computed from calculated fluxes and matched experimental qoz 
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Glucose 



ATP I 



664.3 


NADPH 


142.0 


1.000 




0.214 



i 12.9 i 
i0.6i9i 



i 14.0 1 
i0'02li 



Glucose-6-P 



509.4 



0.767 



6-P-Gluconate 

32.3 



0.049 



NADPH ; C02 



Fructose-6-P 



Ribose-5-P 



ATP 



0.831 

Fructose-1 -6-diP 



Erythrose-4-P 





Gl ycera ldehyde-3-P 

\ 1.118 \ .-^t^^^^ 

1 1.6831 (^^^^^^►nADH 

3-P-Glycerate 



11.1 

017 



ATP 



C02 

73878 



0.209 




1.189 

AcetylCoA 



NADH ; C02 



109.1 
0.164 





112.9 


0.170 


676.8 




1.019 



Oxaloacetate 



NADH 



FADH2 



Succinate 




Citrate 
Isoc^rate 

NADPH ; C02 
Ketoglutarate 




: 29.7 
16.0451 



647.1 



0.974 



NADH ; C02 ; ATP 



; 52.9 
iO.OSO 



q02 
model 


2,210 


3.327 


yC02 
model 


-2,173 


-3.272 



yNADPH 
production 


-961 


-1.446 


qNADPH 

need 


555 


0.836 


yNADH 
production 


-3,343 


-5.033 


qNADH 
need 


0.000 


0.000 



yATP 
production 


-2,657 


-4.000 


qATP 
need 


2,609 


3.927 



mATP 


7,546 


1 (P/0 = 2) 


11.359 



Figure 1. Carbon flux distribution in 5. coelico/or Ml 45. Grey arrows represent carbon fluxes within the central metabolism of S. coelicolor 
IVI145. Arrow widths are proportional to carbon fluxes. Flux values in the central metabolism are included in boxes with solid lines, flux values to 
biomass are included in boxes with dotted lines. The upper numbers represent average actual fluxes (nmol (g dry mass)^' h^'), the lower number 
represent average normalized fluxes (i^mol (nmol glucose)^'). Positive values correspond to consumption and negative values to production. 
doi:1 0.1 371/journal.pone.00841 51 .gOOl 



values very well (difference of less than 2.7±0.5%). Experimental 
and computed CO2 and O2 exchange rates gave respiratory 
quotients {ycm^loz) of 0.982±0.002 for both strains, as would be 
expected during aerobic growth. 

Energy and redox balance 

NADPH production was similar in both strains (Fig. 1, Fig. 2, 
right hand side). However, NADPH requirement for biomass 
biosynthesis was significantly higher for M 1 1 46 while in both, 
NADPH production exceeded needs. The production of NADH 
from catabolism and anabolism was lower (30 ± 3 %) in M 1 1 46 
because the TCA was less active in that strain (see section 
Discussion 4). In Ml 146, the production of ATP from catabolism 
was lower (29 ±2%) than in Ml 45 for the same reason. Contrary 
to Ml 45, the production of ATP from catabolism was insufficient 
to cover ATP needs for biomass synthesis in M 1 1 46 (which are 



13±0.1% higher than in M145, see section 4.3.). Nevertheless, a 
NADH-to-ATP conversion efficiency of 0.25 by the respiratory 
chain in Mil 46 would stiU cover ATP needs in that strain. 

Discussion 

Carbon fluxes in the central metabolism of S. coelicolor 

The aim of this work was to reconstruct a comprehensive model 
of carbon flux distribution of primary carbon metabolism in S. 
coelicolor strains M 1 45 and M 1 1 46, so as to better understand the 
consequences of the absence of antibiotic biosynthesis on central 
carbon metabolism. We used here both Metabolic Flux Analysis 
and ■ C Metabolic Flux Analysis to solve fluxes across central 
carbon metabolism. Our analysis took advantage of published and 
experimental data on composition and reaction stoichiometrics. As 
stated in section Results 4, this model seems valid since observed 
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Glucose 



ATP 



14.6 




0.025 





15.9 
Q.oi.Si 



ATP 



430.7 




ADPH ; CO2 



Fructose-1-6-diP 



0.748 
DHAP 



Glyceraldehyde-3-P 



894.0 





13.0 




0.023 



1.5531 ^ ^^^j^NADH ; ATP 
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Figure 2. Carbon flux distribution in 5. coelicohrM'\'\46. Grey arrows represent carbon fluxes within the central metabolism of S. coelicolor 
Ml 146. Arrow widths are proportional to carbon fluxes. Flux values in the central metabolism are included in boxes with solid lines, flux values to 
biomass are included in boxes with dotted lines. The upper numbers represent average actual fluxes (nmol (g dry mass)^' h^'), the lower number 
represent average normalized fluxes (|xmol (nmol glucose)"^). Positive values correspond to consumption and negative values to production. Boxes 
with a dark grey shading indicate values significantly higher in M1146 than in IVI145. Boxes with a light grey shading indicate values significantly 
lower in IVIl 146 than in M145. 
doi:10.1371/journal.pone.0084151.g002 



respiration rates {qo2 and }'co2 values, Table 4) matched computed 
CO2 release and O2 consumption, with a respiratory quotient of 
nearly 1. 

Nevertheless, Metabolic Flux Analysis is usually viewed as 
unlikely to give a completely faithful representation of the actual 
flux distribution, and consideration of areas in which thecarbon 
flux-balance model diverges from known metabolic behaviour 
may be potentially more informative. For example, the oxidative 
reactions of the oxidative pentose phosphate pathway typically 
carry no flux in flux-balance solutions of heterotrophic metabolism 
in plants [91,92]. This discrepancy arises because mass-balance on 
NADPH can model the provision of NADPH by other possible 
pathways such as, for exemple, isocitrate dehydrogenase. Here, 
isotopic labeling with ['"'C-l] glucose was carried out to infer the 
metabolic commitment to the PPP, which is a source of NADPH. 
There was a very clear ' C-enrichment in TCA-derived amino 
acids [ni] and m2 much larger that they should be at natural ' C- 



abundance), allowing to quantify reliably the ''^C-commitment to 
glycolysis and the TCA versus the PPP. 

Using this approach, we determined a contribution of PPP to 
glucose catabolism in S. coelicolor 1^14:5 (22±0.3%, Fig. 1) close to 
that reported by Naiempoor and Mavituna (25.1%) for S. coelicolor 
during actinohodin production when the growth rate (fi) was near 
0.04 h"' [18]. 

In Ml 146, our estimate of the PPP contribution is 38±0.2% 
(Fig. 2), similar to that found at similar growth rates (ji near 
0.04 h~') in S. lividans (35%) during actinorhodin production [93], 
Streptomyces noursei (33%) during nystatin production [94] as well as 
Streptomjces C5 (30%) during s-rhodomycinone production [39]. 

We postulate that the higher flux through the PPP in Ml 146 is 
the consequence of a higher NADPH demand to sustain both a 
higher growth rate (as previously suggested for S. coelicolor [93,94]) 
and an increased triacylglycerol biosynthesis. Parenthetically, the 
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higher NADPH demand might also contribute to coping with 
oxidative stress (see Discussion 4). 

Reductive power balance 

It is believed that NADPH requirements and production should 
be balanced as cell metabolism reaches the steady state [95] . Here, 
in both strains, NADPH production exceeded NADPH needs 
(Fig. 1, Fig. 2, right hand side). If the NADPH content were 
steady, then the NADPH exces.s would have to be oxidized to 
NADP. In S. coelicolor, this might involve a nicotinamide nucleotide 
transhydrogenase, since our assays showed a thio-NAD reduction 
at the expense of NADPH. In their studies on actinorhodin 
production in an hyper-producing strain of S. lividans, Bruheim 
et al. [96] also suggested the involvement nicotinamide nucleotide 
transhydrogenase to balance excess NADPH production under 
conditions in which actinorhodin was synthesized. 

Similarly, NADH production was larger than NADH needs and 
we assumed that this was compensated for by supplemental 
NADH oxidation through the respiratory chain. Other metabolic 
reactions that oxidise NADH might have also occurred, such as L- 
lactate dehydrogenase or alcohol dehydrogenase. However, such 
enzymatic activities are not likely since bacterial metaboUsm was 
aerobic under our conditions. In addition, we found no L-lactate 
and no ethanol at all in our experiments (data not shown). 

ATP needs, mATP and P/0 

In strain M145, ATP generated by catabolism 
(4.000±0.043 ixmol ATP (jimol glucose)"') would cover ATP 
needs for anaboli-sm (3.927±0.052 |j.mol ATP (|imol glucose)"'). 
To the contrary, an apparent ATP deficit (1.118±0.013 |J,mol 
ATP (|imol glucose)" ) would exist in strain Mil 46. 

For M145 and M1146, 5.643±0.113 jmiol NADH (nmol 
glucose)"' and 4.453±0.089 |imol NADH ((imol glucose)"' were 
generated in excess by metabolism and from putative NADPH 
transhydrogenation, respectively. 

As the carbon balance do not close, we arc- aware that the 
calculation of ATP needs for anabolism may be affected if ATP is 
required for biosynthesis of the unidentified metabolites. This 
should not be the case for NADH, as NADH needs for anabolism 
are negligible (Table 3). 

The amount of ATP generated by NADH reoxidation by the 
respiratory chain obviously depends upon the P/O ratio. Here, 
even a P/O ratio as low as 0.25 would have been sufficient to 
cover ATP needs in Mil 46. 

If P/O ratios larger than 0.25 were taken into account, both the 
oxidation of total NADH to NAD by the respiratory chain and 
intrinsic metabolic phosphorylation (pyruvate kinase, etc.) would 
have produced an apparent ATP excess in both strains. 

This apparent ATP excess, which is certainly balanced by ATP 
consumption for cell maintenance, may be denoted as mATP. 
mATP represents the energy demand for essential, non-growth- 
related processes, including maintenance of ionic potential across 
the cytoplasmic membrane, metabolite transport, cells homeosta- 
sis, macromolecules turn-over and others processes [82]. Since 
80% of the actinorhodin pool in M145 were found in the 
extracellular compartment, a substantial part of mATP was 
probably dedicated to its excretion. 

Calculated mATP strongly depends upon the P/O value used 
and in practice, the P/O value may vary substantially depending 
on growth conditions [97]. However, the P/O ratio chosen here 
has no consequence on the carbon flux distribution found in the 
present study, that is, any value of the P/ O ratio above 0.25 would 
yield the same carbon flux pattern. 



Here, taking into account a P/O value of 2, mATP values (on a 
biomass basis) would be 7.545±0.3 and 4.482±0.2 mmol ATP (g 
cell dry weight)"' h"' in M145 and Ml 146, respectively, that is, 
fell well within the range (4.52-9.64 mmol ATP (g cell dry 
weight) ' h ') of those found in various bacteria [98-102] and S. 
coelicolor [71] grown on glucose as a carbon source. 

A standard way to estimate maintenance energy is to plot 
specific glucose uptake rate (mmol glucose (g cell dry weight) ' 
h ') as a function of the dilution rate Z) (h ') in a chemostat [103], 
the J) intercept being the maintenance energy value [104]. From 
consideration of published dilution rates [10,19,71,93,105], one 
may calculate qs values in Streptomyces and a value of 
0.425±0.05 mmol glucose (g cell dry weight)"' h"' is found. 
Using our own experimental data, maintenance energy should 
correspond to 0.498±0.019 mmol glucose (g cell dry weight) 
h"' inM145 and 0.341 ±0.015 mmol glucose (g cell dry weight)"' 
h"' in Ml 146. Even though such values agree reasonably well 
with the average value from literature, an experimental determi- 
nation of energy maintenance under our conditions remains to be 
determined. 

Alterations of metabolic fluxes in Ml 146 and their 
consequences 

Since Ml 45 produced solely actinorhodin from acetyl-CoA 
during its growth under our conditions, the genetic deletion of the 
cluster for actinorhodin biosynthesis in M 1 1 46 expectedly resulted 
in visible metabohc changes (Fig. 1, Fig. 2). Since acetyl-CoA is the 
main precursor of actinorhodin, acetyl-CoA accumulation could 
have been postulated in Ml 146. Nevertheless, acetyl-CoA 
accumulation is unlikely to occur for a number of biochemical 
[106], physiological [107,108] and thcrmodynamical reasons. 

Other evidences for the absence of acetyl-CoA accumulation 
came from the excretion of larger amounts of pyruvate in Ml 146 
(167±11 |xmol (g cell dry weight)"' h"') than in M145 
(21 ±1 |imol (g cell dry weight) ' h ') during the stationary 
phase. In the absence of growth, acetyl-CoA could not be 
redirected to biomass in Ml 146 and accumulation was simply 
avoided by excreting its direct, freely diffusible precursor pyruvate. 

During the exponential growth of M145, 80.4±1 umol acet)'l- 
CoA (g cell dry weight)" h" are dedicated to biomass produrtion 
and 32.5±1 |lmol acetyl-CoA (g cell dry weight) h to 
actinorhodin production. In Ml 146, with no other modification 
in glucose uptake and carbon flux distribution, a complete 
redirection of 32.5 ±1 |lmol acetyl-CoA to biomass would result 
in a 40% increase of the specific growth rate (0.0552 h '), a 40% 
increase of the anaplerotic flux and a decrease of qo2 (1,581 |J,mol 
(g ceU dry weight)"' h"') and yco2 (1,534 |imol (g ceU dry 
weight) ' h ') values in Ml 146. 

However, the specific growth rate of Ml 146 was only 13 ±0.1% 
higher than that of Ml 45 and the rate of glucose consumption was 
13±0.5% lower. In Ml 146, there was therefore a slightly lower 
flux to acetyl-CoA synthesis due to the lower glucose consumption 
rate and as a result, the flux through the TCA and ATP 
production were lower (Fig. 2). Accordingly, qo2 
(1,475±0.05 (imol (g cell dry weight)"' h"') and ■yco2 
(1,448±0.05 (imol (g cell dry weight)"' h"') values were also 
lower. 

The reasons of the limited increase in M 1 1 46 growth rate are 
not straightforward since S. coelicolor could grow with fi values up to 
0.16 h ', a high value (2 mmol (g ceU dry weight) 'h ') and a 
low actinorhodin production rate {yAcr— 6 nmol (g cell dry 
weight)" h" ) [19]. Three mechanisms at the origin of growth 
restriction are plausible: 
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First, in Streptomyces peucetkus var. caesius and S. coelicolor, a 
transcriptional regulation of glucose uptake rate has been 
demonstrated [36,109,110]. In Ml 146, a limitation of glucose 

entry mediated by an as yet unidentified regulatory protein (with 
acetyl-CoA or triacylglycerol as co-regulators) might have 
occurred. 

Second, the absence of actinorhodin synthesis in Ml 146 might 
also cause growth limitation as a side effect. Indeed, quite recently. 
Shin et al. [Ill] demonstrated that the extracellular pool of 
actinorhodin was involved in the activation of the SoxR regulon in 
S. coelicolor. This regulon encodes proteins involved in the transport 
or the turn-over of small redox-active toxic compounds [112]. As 
described in E. coli by Fridovitch [113], such compounds might 
alter catalytic [4Fe-4S] centers of various dehydratases such as 
fumarase and aconitase in the TCA, or the dihydroxyacid 
dehydratase involved in the synthesis of abundant amino acids 
(lie, Leu and Val), pantothenic acid and coenzyme A. Thus, under 
70% of dissolved oxygen, the lack of actinorhodin production in 
M 1 1 46 might have lead to an increased intracellular oxidative 
stress, which in turn impeded the TCA, ATP synthesis, biomass 
synthesis and glucose consumption. The involvement of oxidative 
stress would be consistent with the increase in NADPH production 
by the PPP in Ml 146 (see Discussion 2). 

Third, the carbon metabolism in Ml 146 pardy compensated for 
the reduced growth rate (lower /j) and antibiotic synthesis by 
a[:[:umulating neutral lipidic compounds (triacylglycerol, see 
Results 1). In a similar molecular context, that is, in the bldA 
mutant of S. coelicolor unable to produce actinorhodin and 
undecylprodigiosin, Plaskitt and Chater [1 14] also demonstrated 
a higher triacylglycerol content. In fact, there seems to be a 
compromise between antibiotic and triacyglycerol production 
since (i) a higher actinorhodin production has been observed in 
the sco0958 mutant impaired in triacylglycerol synthesis [115] and 
(ii) the simultaneous deletion of the gene clusters encoding 5- 
polyketide synthase and addition of a plasmid comprising the gene 
cluster encoding actinorhodin biosynthesis, resulted in a dramatic 
increase of actinorhodin production in S. coelicolor [116]. 

Ml 146 as a super-host 

Ml 146 was successfully used as a super-host for the production 
of actinorhodin, congocidin, chloramphenicol and caprazamycin 
after introduction of the corresponding gene clusters [10,11]. 
However, the use of Ml 146 did not allow the efficient production 
of clorobiocin and coumermycin Al [10]. 

The metabolic pathways associated with congocidin [117], 
chloramphenicol [118], caprazamycin [119], clorobiocin and 
coumermycin Al production [120,121] have all been elucidated. 
Precursor metabolites for the production of congocidin (3- 
phospho-glycerate, N-acetyl-glucosamine, sedoheptulose- 7 -phos- 
phate and cytosine), chloramphenicol (shikimic acid, dichloro- 
acetyl-CoA), and caprazamycin (uridine, glycine, S-adenosyl- 
methionine, 3-methyl glutaryl-CoA, fatty acyl-CoA, per-methyl- 
ated rhamnose) are all derived from glucose- 1 -phosphate, mole- 
cules from the pentose phosphate pathway, 3-phospho-glycerate 
and acetyl-CoA. 

In the case of clorobiocin and coumermycin Al, precursors are, 
in that order of incorporation, prephenate (derived from both 
erythrose-4-phosphate and phospho«no/pyruvate), noviose (derived 
from glucose- 1 -phosphate) and 3- or 5-methyl-pyrrole-2-carbox- 
yhc acids. These pyrrole groups originate from threonine 
(precursor oxaloacetate) or proline (precursor c/.-ketoglutarate). 
Flinspach et al. [10] have clearly shown that clorobiocin and 
coumermycin Al production was unsuccessful! in Ml 146, the 



main products beeing aminocoumarin derivatives without any 
pyrrole group or not the expected ones. 

Although we are aware of possible pleiotropic effects due to the 
introduction of f()r(;ign genes in a given genetic background, these 
results seem consistent with our carbon flux description in Ml 146, 
i.e., an enhanced flux through PPP, a moderately decreased 
glycolytic flux to acetyl-CoA, an increased acetyl-CoA extraction 
for triacylglycerol synthesis and a lower flux in the TCA. In 
clorobiocin and coum(;rmycin Al biosynthescs, the consumption 
of er\''throse-4-phosphat(' and pliosphoajo/pyruvate may exacer- 
bate the flux impairement through the TCA, leading to the 
impossibility to produce a sufficient pyrrole content. 

Conclusion and perspectives 

In this paper, we demonstrated that combining Metabolic Flux 
Analysis and ' ^C Metabolic Flux Analysis yielded an accurate 
description of carbon fluxes in S. coelicolor under steady-state 
conditions. Our work further suggested the key role of nicotinamide 
nucleotide transhydrogenase for redox homeostasis in S. coelicolor. 

The main metabolic difference between Ml 45 and Ml 146 
consisted in an enhanced flux through the PPP in Ml 146, likely to 
face higher needs for NADPH in this strain. The metabolic 
differences between Ml 45 and Ml 146 also highlighted the 
energetic cost for actinorhodin production and excretion in Ml 45. 

In other words, there would be a metabolic trade-off between 
instantaneous growth and the production of antibacterial com- 
pounds, which are likely to provide a competitive advantage [122] 
or to be involved in quorum sensing [112,123,124]. In other 
words, antibiotics are beUeved to help to defend the food source 
and to face bacterial competition in natural environments. 

Our results also suggest a plausible competition between 
polyketide biosynthetic pathways for their common precursor, 
acetyl-CoA. Therefore, the deletion of polyketide antibiotic 
pathways might be instrumental in increasing acetyl-CoA avail- 
ability and thus triacylglycerol, which is an important source of 
renewable oleochemicals. Conversely, the impairement of triacyl- 
glycerol synthesis might be of interest to improve polyketide 
antibiotics production yield. 

In Streptomyces, antibiotics are usually produced in small amounts 
at the transition phase when bacterial growth slows down as a 
result of nutrient exhaustion, that is, not in steady-state conditions. 
New rapid sampling techniques [125], new modelling approaches 
and highly sensitive mass spectrometers make possible to carry out 
non-steady '''C Metabohc Flux Analysis [126,127]. Together with 
our knowledge of S. coelicolor metabolism under steady-state 
conditions, it is now possible to study metabolic carbon fluxes 
during non-exponential growth stages. This will offer a better 
insight into the understanding of secondary metabolite biosynthe- 
scs under such conditions. 
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